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Abstract ' • 

i • • 

The structure of the drift instability is examined in several den- 
sity regimes. Let X be the total electron mean free path, k the wave- 

G Z 

vector component along the magnetic field, and v x /v|| the ratio of perpen- 
dicular ion diffusion to parallel electron streaming rates. At low densities 
(k > 1) the drift mode is isothermal and should be treated kineticly. 

In the finite heat conduction regime ,/m/M" <k X < 1 the drift insta- 

z e 

bility threshold is reduced at low densities (v_l/v|| < 0.1) and increased 

at high densities (v x /v|| > 0,1) as compared to the isothermal threshold. 

Finally, in the energy transfer limit (k X < / m/M ) the drift; instability 

z e 

behaves adiabatically in a fully ionized plasma and isothermally in a par- 
tially ionized plasma for an ion-neutral to Coulomb collision frequency ratio 

V . /v . . > 2 .Am / M at T . = T - T . 
in u i e n 
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Introduction 

i 

The most obvious source of free energy for the instability of the 
nighttime equatorial F region called spread F is the sharp plasma density 
gradient on the underside of the postsunset F layer. While a variety of 
density gradient dependent modes have been considered (cf. Farley et al > 

1970) to explain the geomagnetic field alligned density irregularities with 
perpendicular wavelengths from a few m to km, Hudson and Kennel (1974a) 
are the first to consider collisional drift waves driven solely by the density- 
gradient. The collisional drift mode appeared in the electrostatic limit of the 
electromagnetic dispersion relation for low frequency density-gradient driven 
modes in a partially ionized collisional plasma derived in the companion paper 
(Hudson and Kennel 1974b). The observed altitude threshold for spread F on- 
set (Farley et al . 1970) suggested the extension of previous work on the col- 
lisional drift mode (Tsai et al . 1970) to partially ionized plasmas where 
neutral collisions affect the instability threshold and growth rate. 

Recent work on the collisional drift mode has shown that its growth 

rate is modified by finite electron heat conduction along the magnetic 

field at parallel wavelengths exceeding the electron momentum transfer 

mean free path \ . However, there is an apparent paradox in the litera- 
c 

ture; Ellis and Motley (1974) find that the drift mode growth rate is enhanced 
by finite heat conduction, while Tsai et al . (1970) conclude that it is reduced, 
especially at high densities. We resolve this conflict by obtaining an ex- 
plicit criterion which separates the finite heat conduction regime /rn7 M < 
k X <1 into low and high density regions where the drift mode growth 

ZG 
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rate is enhanced or reduced from its isothermal value. A single parameter 
which is the ratio of the perpendicular ion diffusion to parallel electron 
streaming rates, v x /v| f - 0.1, separates the low and high density regions. 

In application to equatorial F region irregularities v / v n < 0. 1 in general (Hud 
son and Kennel 1974a) and the drift mode growth rate is enhanced by finite 
heat conduction. 

% ..... 

Since equatorial F region geometry permits very long parallel 

wavelengths, we have extended the analysis of the collisional drift mode to 
parallel wavelengths exceeding the energy transfer mean free path, k X < 
/rn/M , where we find that temperature fluctuations and finite heat conduc- 
tion effects are surpressed. In this long wavelength regime which we call 
the energy transfer limit, a fully ionized plasma behaves adiabatically, 
because there is a balance of energy exchange between elecUons and ions. 

A partially ionized plasma with v. /v.;> 2 /"in JM at T - T. - T behaves 

m n e i n 

isothermally (Hudson 1974), because energy in electron temperature fluc- 
tuations is lost to the ions and then lost to the neutral sink faster than it 
can be transferred back to the electrons. 

Electrostatic Dispersion Relation 

Finite electron heat conduction, energy transfer between species 
and neutral collisions have all been included in our derivation of the elec- 
trostatic dispersion relation for the collisional drift mode from two fluid 
theory in the companion paper (Hudson and Kennel 1974b) 


4 . 
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neglecting gravity which has a minor effect on the drift mode and acoustic . 

terms which are stabilizing for k C > (Hudson 1974). Here L ± is 

z s D 

2 2 2 

the density gradient scale length, b = k C /n is the ion FLR parameter, 

y s 

m --/BC -/L* is the diamagnetic drift frequency in terms of the ion 
I ) s 

2 

acoustic speed C g = T^/M, vj_/b = v + 0. 3 b is the perpendicular ion 

2 2 . 

diffusion rate, vj| = is the parallel electron streaming rate. 


v - v . + v _ is the total electron collision frequency, v.. and v. (j = 3, i) 

' / T , 

T * i e 


c ei es 

arc Coulomb and neutral collision frequencies . p = 1 + R = 1 + T./T . 


Vve also define . . 

X - 2/3tC r C x *(l*C t ) ! ]*l!(m.«)» e ^, • 

' - «j(. ♦£,)* .s(. « t )(. - « t ) , , . (, . R T ) w„|( Ven - 3v.,)/v, 

which are combinations of finite electron heat conduction terms, plotted 

in Figure 1 neglecting the energy transfer terms proportional to m/M 

which dominate at high densities and long parallel wavelengths (k \ < /m/M ). 

z e 

In the electrostatic dispersion relation (la) the drift mode is coupled 
to the entropy mode (Tsai, et al . 1970; Hudson and Kennel 1974b), interchange 
and acoustic terms. Finite parallel electron heat conduction is included 
in (la), which reduces to the electrostatic dispersion relation obtained by 
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Tsai et al . (1070), when neutral collisions and energy transfer from elec- 
trons to ions and neutrals are neglected. 

Isothermal Limit 

The isothermal limit which applies to k X >1 is formally recovered 

Z 6 

from (la) by letting the dimensionless coefficient of electron thermal con- 
ductivity C -* ®. Hence x -* 03 in contrast to § which remains finite, 

!x» 

and (la) yields 

b{( a + R-j-Wq) (w + I VjL /b ) + g/L A } 

* * 1 v»fl C«C 1 + pb] - U) D + Iv A p) (lb) 

The collisicnlc cz limit k \ •> 1 is only fully treated by a kinetic theory 

(Coroniti and Kennel 1970) which includes Landau resonance effects. 

Neglecting these resonances (lb) reduces to the fluid limit of their kinetic 

dispersion relation as 0 or V|| 

In the opposite limit v n -♦ 0 for H r ^ ~ 1, we again find that X -* « 

while £ remains finite. Hence by including energy transfer effects in X 

and 5, we have obtained the short parallel wavelength isothermal results 

at long parallel wavelengths. Including ion temperature fluctuations neglected 

in the derivation of (la), we have found that. this result is valid in a partially 

ionized plasma with v. /v.. > 2 / ml M at T * T. = T (Hudson 1974). 

in n e i n 

For this small ratio of ion-neutral to Coulomb collision frequencies, the 
ion energy loss rate to the neutral sink exceeds the ion energy transfer 
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rate back to the electrons, and the plasma behaves isothermally at parallel 
wavelengths exceeding the electron energy transfer mean free path /M/m \ « 
However, in a fully ionized plasma there is a balance of energy exchange . 
between ions and electrons, and the plasma behaves adiabatically at long 
parallel wavelengths. 

The isothermal dispersion relation (la) has therefore acquired a 
new significance at long parallel wavelengths where Landau resonance 
effects are negligible. Note that the cubic dispersion relation (la) reduces 
to a quadratic (lb) with the elimination of the entropy mode in, the isother- 
mal approximation. It is useful first to obtain the marginal stability criterion 
in this isothermal approximation for the simple quadratic, since the same ' 
proceedure will be used to derive the general marginal stability criterion 
for the cubic. ’ 

Isothermal Marginal Stability Criterion 


At marginal stability w is real and the real and imaginary parts 
of (lb) vanish separately. The imaginary part yields the dispersion relation 


(l - R T Vj/v«) 

® * W D (UpbVv A 7^7T* 


1 + R. 


(2a) 


For H^Vx/vji < < 1 we obtain the drift mode dispersion relation 


u) r: uj^/d + pb) 


(2b) 
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The drift mode marginal stability criterion is obtained by substituting 
(2a) back into the real part of (lb), neglecting gravity 


* ( 1 ~ v») 

(1 + pb + vj_/ \>\\ ) 




. 2 . 2 

*k z L x 


(3) 


Instability results for 0 < K < K, where K is the destabilizing parameter 
proportional to the density gradient scale length Lj.. An unstable drift 
mode is only obtained for R^vx/vjj < 1, which requires finite k^. Equation 
(3) is plotted in Figure 2 for R T - 1 and b < < v x /v„ < 1. The threshold 

a 2 2 

increases (Ka k z L x ~ decreases) with increasing Vj_/v|j proportional to 
total density. 

General Marginal Stability Criterion 

Having analyzed the marginally stable quadratic (lb) containing 
the FLR coupled interchange and drift modes in the isothermal limit, we 
now consider the marginally stable cubic (la) containing the FLR coupled 
interchange, entropy and drift tpodes at intermediate parallel wavelengths 
/rn/M < k X <1 where finite heat conduction effects apply in general. 

The real and imaginary parts of (la) vanish separately at marginal 
stability where y) is real. The imaginary part, which was linear in the 
isothermal case, is now quadratic with the addition of the entropy mode. 
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l Np.imiim; that quadratic assuming | 1 - v x /v ,| j > b (3R T + K) yields 
{1n> drift mode root 

m + ( 1 ^ R T v -t / v 0 " ^TX b W ~ 5 )(l + R |)bK 

‘"O [ 1 + Vj/vh + ( 1 + Rj + x) b ] (l-R T v x / v n) 

(4) 


The instability criterion for the drift mode is obtained by substi- 
tuting (4) back into the real part of (la), yielding (Tsai etal., 1970) 



- x<vj,/vx><x ■ 5) + *{Me - o(i +r t ) + fj 



(5) 


This is the most general drift mode marginal stability criterion which 

applies to all densities k \ <1, 

z e 

We have plotted the marginally stable K vs. v x /v,| in Figure 3 for 
a Jully ionized plasma at « 1 and several values of the parameter N = 0. 3 
proportional to density. We use the fully ionized values of X - 3. 03 
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and ?, ~ 1.95 (Tsai et al . 1970) and substitute for b = (l/N J Kv^/vjj . Com- 
parison with the isothermal curve in Figure 2 demonstrates that the drift 
instability threshold is lower at low densities (v^/v^ < < 0. 1) than in the 
isothermal limit. 

Analytic approximations to the high and low density limits can be 
obtained from the left and right hand sides of (5) respectively. At high 
densities the left hand side of (5) yields (Tsai et al . 197J)) 


( 6 ) 


which is plotted as the horizontal dashed curve in Figure 3. At low den- 
sities the right hand side of (5) yields 




(8a) 


which approaches 

K<X/(X -?> (8b) 

as Vj./vu - 0. Equation (7) is plotted as the vertical dashed curve in 
Figure 3. These approximate marginal stability criteria again demonstrate 
the increase at high densities and decrease at low densities, as compared 
to the isothermal threshold (3). 



10 . 


It is significant to note in Figure 3 that for all Tsl values the exact 
marginal stability curves (5) cross at v x /v„ = 0. 1, K = 0.25 where the 
approximate curves (6) and (7) cross. This single parameter point 
separates the finite heat conduction behavior of the drift mode into high and 
low density regions where the threshold is increased or decreased from 
its isothermal value. 

, In both the large and small v x /v|| limits the finite heat conduction 

marginal stability criteria (6) and (7) reduce to the isothermal threshold 

(3) as X In (6 ) (1 + £ ) comes from 1 + y - (y - £ ), and 5 remains 

finite as we take the isothermal limit’x -* 00 . Hence for k \ < / rn / M the 

z e 

plasma behaves isothermally, independent of v 1 /v t!# , under the assumptions 
discussed in the preceeding section. 

Drift Mode Growth Rate 

Assuming u>/v,|, v ± /v,j and bR^^/v,, < < l decouples the drift mode 
from entropy and interchange mode terms in (la) and yields an approximate 
expression for the drift mode dispersion relation and growth rate in the 
low density finite heat conduction limit 

a = <d d /(1 +pb) (15a) 

Y/ V i - p/( 1 + pb )[“/ ( x - f)(K' '/(l + 0 b)j - 1 j 


(15b) 
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f > 

Neglecting FLR corrections, this result agrees with Ellis and Motley (1974). 
In the limit X -* » the isothermal growth rate is recovered (Chu et al. 1969). 
It is apparent that finite heat conduction increases the drift mode growth 
rate by the factor X/(X - 5) = (5/3) in a fully ionized plasma at low densi- 
ties (v j_ / v o' < 0. 1). At higher densities(v 1 /v) ( > 0. 1) and well above marginal 
stability where the preceeding assumptions break down the drift mode growth 
rate must be obtained by solving the cubic dispersion relation (la) in general. 

The drift mode propagates in the electron diamagnetic drift direction. 
It is destabilized by the difference between the electron guiding center 
current and the ion current which is an ion FLR effect of order b. In con- 
trast, in the companion paper we have seen that the interchange mode is 
stabilized by FLR effects in a weakly collisional plasma (v /b < 2/g/L x ). 

We therefore expect the drift mode growth rate to peak at shorter perpen- 
dicular wavelengths (Hudson and Kennel 1974c). Damping is proportional 
to the perpendicular ion diffusion rate Vi/b - + 0.3 by... 

Conclusions 

We have obtained a marginal stability criterion for the drift mode 

in a partially ionized plasma, which is valid for all wavelengths greater 

than the total electron mean free path (k \ < 1). An apparent contradiction 

z e 

in the effect of finite heat conduction on the instability threshold and growth 
rate has been resolved (Tsai et al . 1970; Ellis and Motley,. 1974). At low 
densities {v x /vj| < 0.1) the drift mode decouples from the entropy mode 
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and the growth rate of the former is enhanced. At higher densities 

(vx / v if > 0. 1) the drift and entropy modes are coupled, and the drift mode 

growth rate is suppressed. 

When k X < fmj M the drift mode approaches the isothermal 
z e 

electron limit for all v ± /v||, and we recover the isothermal instability 

criterion of Chu et al . (1969) generalized to a partially ionized plasma. 

This assumes that ions behave isothermally due to energy loss to the neutral 

sink, which Hudson (1974) has shown is valid for v. /v.. > 2 /m/M at 

in n 

T = T = T In a fully ionized plasma there is a balance of energy ex - 
e i n' 

change between ions and electrons in long parallel wavelength modes 

(k X < J m/M) which behave adiabatically . This distinction in behavior 
z e 

at long parallel wavelengths is significant since the drift instability threshold 
and growth rate is modified by the nonisothermal effects considered above. 
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Figure Captions 

Figure 1. Plot of energy heat flow parameters for R^ = 1 neglecting 

energy transfer terms. X = 2/3 [C C, + (l + C ) 2 ] and § -2/3 (1 + C ) 2 , 

r x t t 

where C , C., and C are the dimensionless resistive, thermal force, 
r t x 

and thermal conductivity coefficients computed by Braginskii (1965) for a 

fully ionized plasma: C = 0.15, C .= 0.71 and C * 3.16. Schunk 

r t x 

and Walker (1970) have plotted them for a partially ionized plasma as 

g = C , g fg = C, 1 and (l/g, -0,5g /g g ) = C / 5 from Shkarofsl<y 's 
rr r cr ti. t k cur x 

(1961) tabulation of the g-factors. Here we assume that electrons collide 
with hard sphere neutrals at a rate proportional to the electron thermal 
speed v = a e . 

Figure 2. Plot of the isochermal drift mode marginally staple parameter 
2 

K = v x v]| /bu>^ as a function of v x /vn from (3) at s 1 for b < < v x /v|, < 1. 

a 2 2 

The threshold increases (K « k L x decreases) with increasing v x /v|j 

z 

proportional to total density. 

Figure 3. Plot of the fully ionized finite heat conduction drift mode 

2 

marginally stable parameter K s v x v|| /bou^ as a function of v x A>[) from 

(5) at = 1 for several values of N = 0. 3 L x /x e proportional to density. 

The exact curves are well approximated by the dashed lines (6) and (7 ) 

at high and low densities; they correspond to the left and right hand sides 

of (5) which are separately independent of N. The crossover point at 

v x /v(] ~ 0.1 corresponds to the separation between the low and high density 

finite heat conduction regimes where the drift instability threshold is lowered 
A 2 2 

(K «:k Lj_ increased), and raised compared to the isothermal threshold in 

■ / , 

Figure 2. ' 
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